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■r::j- ' Models of strongly correlated electrons that tend to phase separate are 

studied including a long-range 1/r repulsive interaction. It is observed that 
a \ 

Vh ' charge-density-wave states become stable as the strength of the 1/r term, 

P^ ' Vcoui, is increased. Due to this effect, the domain of stability of the supercon- 

o' 

J-^ \ ducting phases that appear near phase separation at Vcoui = is not enlarged 

L^ by a 1/r interaction as naively expected. Nevertheless, superconductivity ex- 

d ' ists in a wide region of parameter space, even if phase separation is suppressed. 

Our results have implications for some theories of the cuprates. 
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The presence of charge modulations and phase separation in the high-Tc compounds, as 
well as in other related materials, has recently attracted considerable attention. Experimen- 
tal evidence has shown that the 38K superconductor La2Cu04+y (LCO) has a regime where 
phase separation occurs, 0] which is adjacent to a superconducting phase. Hole-rich and 
hole-poor phases coexist at temperatures below ~ 230K. While a similar macroscopic phase 
separation does not occur in La2_xSrxCu04 (LSCO), a microscopic segregation of doped 
holes forming walls of hole-rich material separated by undoped domains has been reported. 
Phase separation has also been observed in nonmetallic, nonsuperconducting La2Ni04+y 
(LNO), while electron diffraction techniques applied to La2-xSrxNi04+y have indicated 
the presence of charge modulations at low temperatures. ||^ On the theoretical side, the 
description of electrons and holes in the high-Tc cuprates is still under much discussion. 
Among the many proposed theories, scenarios linking phase separation with superconduc- 
tivity have been discussed. [^^ In the widely studied two dimensional t-J model phase 
separation exists at large J/t, and in its vicinity strong superconducting correlations have 
been detected in the ground state. |^0 Actually, the presence of superconductivity close 
to phase separation in electronic models seems a general phenomenon. [^ Fig.l(a,b,c) illus- 
trates these ideas for the particular cases of the one-dimensional (ID) and two-dimensional 
(2D) t-J model, and also for the 2D t-U-V model i.e. the standard Hubbard model with 
a nearest-neighbor density-density interaction. There is a growing perception that phase 
separation and superconductivity are intimately related in electronic models. PJ^ 

While a scenario linking phase separation and superconductivity is appealing, the accu- 
mulation of charge in the hole-rich phase makes relevant the long-range repulsive Coulomb 
interaction among the carriers. Certainly such interactions will render the phase separated 
regimes unstable in Fig.l, while the superconducting phase, being uniform, should not be 
strongly affected. Thus, a natural interesting possibility arises: can it occur that once 
phase separation (PS) is suppressed by a long-range force, the parameter regime formerly 
occupied by PS in Fig.l(a,b,c) is taken over by superconductivity? This effect would no- 
tably enhance the region of stability of the superconducting phase. Actually, Kivelson and 



Emery predicted that high temperature superconductivity is estabhshed once long-range 
Coulomb interactions eliminate phase separation. However, preliminary studies in the 
ID t-J model monitoring the pairing correlation functions showed that superconductivity 
was not enhanced by the 1/r term. |^ It was there conjectured that charge-density-wave 
(CDW) phases may be responsible for this effect. More recently, a study of a classical spin- 
one lattice gas model including long-range repulsive interactions among spins observed that 
the equivalent of ordered charge-modulated states become stable. |T^ Similar indications of 



CDW phases have also been reported in the context of the three band Cu-0 chain. [|TI 

The purpose of the present paper is to discuss a detailed analysis of the competition 
between superconductivity, CDW and phase separation in electronic models including 1/r 
repulsive forces. We considerably improve upon previous calculations by using realistic 
models with mobile electrons not only in ID but also in 2D, and explicitly monitoring 
the correlation functions corresponding to the three competing phases. More specifically, 
the t-J model have been studied using exact diagonalization techniques, and the 2D t-U- 
V model, an effective model for d-wave superconductors, has been also analyzed within 
a mean-field approximation. We conclude that even in the presence of mobile carriers, 
CDW phases indeed are stabilized by the 1/r term in a large region of parameter space, 
reducing the potential domain of stability of the superconducting phase. Nevertheless, we 
have observed that superconductivity is not strongly suppressed either (unless Vcoui crosses 
some threshold), and thus the fluctuating phase separation scenario for the cuprates, with 
some modifications to account for the competition with CDW, remains a viable possibility. 
The strong competition between superconductivity, CDW, and phase separation in electronic 
models seems a general phenomenon. 

The Hamiltonian for the t — J model with 1/r interactions has the form 

Htj = -t 51 (CivCja + h.c.) + J 51 (Si • Sj - -UiUj) 

<iJ>(T <i,j> 

ij J^ij 



where the c-operators are hole operators acting on non-doubly occupied states, ry is the 
shortest distance between sites i and j, and the rest of the notation is standard. To analyze 
the ground state properties of this Hamiltonian exact diagonalization techniques are used. 



|T^ Charge, spin and pairing correlations are monitored as a function of J/t, Vcoui/t, and 



(n). In ID, the conformal field theory parameter kp is also studied. |TB|. This parameter 
is defined as kp = 7fVc{n)'^K,/2, where the charge velocity Vc is obtained from the energy of 
the lowest spin-singlet state in the subspace with momentum k = 2tt/N (N is the number 
of sites), while the compressibility k at a given number of particles Ng is obtained from the 
ground state energies with Ne, Ne + 2 and Ne — 2 electrons, following standard procedures. 
T3|,|9[| li kp > 1 the singlet pairing correlations decay the slowest against distance, and 



thus are dominant in the ground state. To discuss the results in both the ID and 2D t-J 
model, for simplicity we will use the quarter filling density, (n) = 1/2, where the signals 
of superconductivity are clear even in 2D small clusters. We have checked that for other 
densities (where superconductivity exists near phase separation at Vcoui = 0), the results are 
qualitatively similar as at (n) = 1/2 (but with CDW states having a different modulation). 
To start the analysis, it is helpful to consider the atomic limit Vcouh J 3> t, where intu- 
ition can be gained about the states that will compete with superconductivity and phase 
separation. This limit can be studied accurately on large clusters and as anticipated from 
the spin one model, 0] several CDW phases were observed in the ground state as a function 
of J/Vcoui after the 1/r interaction makes the phase separated regime unstable. These CDW 
phases have an increasing number of electrons in each microscopic cluster as J/Vcoui grows, 
since J favors the formation of large spin structures to gain antiferromagnetic energy (see 
Fig. 2). Monitoring the density-density correlation functions, we observed that the CDW 
states are stable even for a finite hopping t, and their rough domain of stability is shown 
in Fig. 2. Phase I is a standard Wigner crystal. Phase II is a Wigner crystal of pairs i.e. a 
regular distribution of charge 2e spin- singlets, similar to those observed in the t-J-V model. 
I^J^ This state is stable since the pairs take advantage of the short range effective attractive 
force produced by J. Phase III has clusters with four electrons, and as J increases the size 



of these microscopic clusters also increases smoothly producing a cascade of CDW phases. 
In the limit where J is the only relevant scale, phase separation is recovered. |]1^ 

When the hopping amplitude t is nonzero, intermediate values of J, and a small coupling 
Vcoui, a regime of superconductivity exists in the phase diagram as shown in Fig. 2. Its 
boundaries were obtained using exact diagonalization techniques on chains of 12 and 16 



sites. ||T2| This study shows that superconductivity dominates (i.e. kp > 1) on a finite region 
of parameter space, while in the bulk limit phase separation exists only along the Vcoui = 
line in Fig. 2. In the pure t-J model, the compressibility diverges at the boundary of phase 
separation, and thus kp on a finite cluster changes sign as shown in Fig. 3a. On the other 
hand, at finite Vcoub ^p is a smooth function of J/t, it becomes larger than one on a small 
region, and then it smoothly converges to zero at large J/t. To gain more intuition about 
the physical behavior of the system, we also studied pairing correlations observing that in 
the regime where kp > 1, these correlations are indeed very strong in the ground state and 
they continue having a large value, even up to Vcoui/t ~ 0.5 (with J/t ~ 3) i.e. well beyond 
the apparent stability regime signaled by kp. This curious effect shows that short distance 
superconducting fluctuations may be relevant in a wide region of parameter space, even if 
their asymptotic power-law decay is not the dominant one. On Fig. 3b, the superconducting 
(pairing) and CDW (charge) correlations (phase IV of Fig. 2) summed over all distances are 
plotted at a flxed Vcoui = O.lt against the exchange coupling. The competition between the 
two phases is clear showing explicitly that even when phase separation is made unstable, 
the region left empty is not necessarily taken over by superconductivity but by a crystal-like 
charge modulation. Thus, our results conflrm the potential importance of charge-density- 
wave phases in electronic models for the cuprates. ||T0|,|9|,pJ| 

The analysis of the 2D t-J model is more complicated since there is no analog of kp, and 
the linear cluster sizes are smaller than in ID. Nevertheless, 2D is the dimension of interest 
for the cuprates, and it deserves to be explored. We have carried out an analysis of this 
case using 4x4 clusters, and monitoring the pairing and charge correlations to set up the 
phase diagram. In the atomic limit, which can be explored on larger clusters, the set of 



stable CDW phases is qualitatively different from the ID case. As observed before in the 



context of a classical spin-one system |T^ and for the t-J model close to half-filling |]T5 
striped phases (holes ordered in one dimensional chains along the x or y axis) are dominant 
in most of parameter space. As J increases, the number of contiguous chains of electrons in 
the striped phases smoothly increases. Away from the atomic limit, a numerical analysis of 
the pairing correlations similar to that carried out for the pure 2D t-J model (see Ref. for 
technical details) suggests that superconductivity is robust in a region of parameter space 
analogous in shape to that observed in ID. While the actual boundaries between the various 
phases in Fig. 3c should be considered only as qualitative, nevertheless the trends are clear 
from the numerics, and similar to those in ID. Thus, our study suggests that although phase 
separation is destabilized by the 1/r correlation, the region left empty is not only replaced 
by the neighboring superconductivity but also by CDW phases. 

To show that this effect goes beyond the t-J model, and it may actually affect all elec- 
tronic models with coexisting superconductivity and phase separation at Vcoui = 0, here we 
consider the t-U-V model. Its Hamiltonian is 

H = -t E (Ci^Cj. + h.c.) + Uj^Ui^n;^ +VY. ni'^J 

<i,j>o- i <ij> 

+VeoulE— ' (2) 

ij ^ij 

where the notation is standard (for details see Ref. [Q). The phase diagram at Vcoui = 
and half-filling was shown in Fig.lc using a mean-field approximation. For simplicity, 
let us consider U/t = which is at the center of the stability region of the d-wave su- 
perconducting phase, and an attractive V term which can be effectively produced by 
spin-polaronic correlations in a more realistic Hamiltonian (results at other values of 
U/t are qualitatively similar). As was shown before, |^ this model at half-filling has 
a stable d-wave superconducting phase with order parameter Aq = -^(cosqx — cosqy), 
where Aq is obtained through the self-consistent equation X^k /^^ 2 ~ ~ ^P' "with 
Veflf(q) = 2^ Er exp(iqr)[V(5r,x + 5r -ic + (^r,y + ^r -y) +^^]. The grouud state energy of 



this condensate is Eqs = ^/y^ — y " ^ I^k v ^k + ^k + %f 5 where Eps is the energy of the 

6 



phase separated state (half the lattice empty, the other half doubly occupied). The en- 
ergies of the CDW and phase separated states can also be written in closed form, but it is 
impractical to quote them here. We have also introduced a cutoff of four lattice spacings 
for technical reasons dealing with the 1/r interactions. This makes phase separation stable 
on a finite region for any finite cutoff. A straightforward comparison of the energies of the 
competing states was enough to establish the mean-field phase diagram shown in Fig.4a. As 
it happened for the t-J model, in this case also the d-wave phase is eventually suppressed 
by CDW states. Note that there is a large region of stability for the d-wave phase but 
this may be a consequence of the approximations made in evaluating the energy of the dif- 
ferent phases. However, we do not expect qualitative changes in the phase diagram given 
improved approximations. Note that the effect of CDW and phase separation is usually 
neglected in studies of models of d-wave superconductors, which Fig. 4b suggests to be a 
dangerous practice. 

In this paper, we have studied the competition between superconductivity, phase sep- 
aration and charge ordered states in models of correlated electrons both in ID and 2D. 
Increasing the strength of the 1/r repulsive interaction, we found that superconductivity 
does not take over the entire region left empty by phase separation, since CDW states be- 



come competitive in energy as observed in the spin-one lattice gas model. ||IOl Nevertheless 
the superconducting phase survives in the presence of the Coulomb interaction as long as its 
strength is not too large. In Fig. 4b, we summarize our results with a schematic representa- 
tion of the influence of long-range interactions over an idealized system which at Vcoui = 
has both superconductivity and phase separation. "g"is the strength of a generic coupling 
constant which produces an attractive short-range force in the system. It is interesting 
to attempt to locate some of the known cuprates in different places of this phase diagram 
(Fig. 4b) . LCO seems to have a small effective Coulomb interaction or screening length since 
phase separation exists in this compound, while LSCO may have intermediate values of 
these parameters that suppress phase separation. LNO (due to its small mobility caused by 
the large Ni spin) seems located in the CDW region of Fig. 4b. Summarizing, in electronic 



models where effective sliort-range attractive forces are operative, superconductivity com- 
petes for stability with phase separation if the repulsive 1/r interaction is neglected, and 
with CDW states if 1/r interactions are included. 
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Figure Captions 

1. Illustration of the presence of phase separation near superconductivity in electronic 
models without 1/r interactions: (a) Phase diagram of the ID t-J model as a function of 



the electronic density {n) and J/t from Ref. [|T3]|. PS denotes "phase separation", kp> 1 



{kp < 1) implies that singlet superconducting (spin and charge) correlations decay the 
slowest in the ground state; (b) Schematic phase diagram of the 2D t-J model taken 
Ref. [§]. The notation is similar to (a) with x = l — {n), "SC" denoting a mostly d-wave 
superconducting region, "PM" a paramagnetic region, and "AF" antiferromagnetism; 
(c) Mean-field phase diagram of the 2D t-U-V model at half-filling (Ref. §]). SS, 
DS, CDW, SDW, and PS denote s-wave superconductivity, d-wave superconductivity, 
charge- density- wave, spin-density-wave, and phase separation, respectively. 

2. Phase diagram of the ID t-J model with long-range 1/r interactions (see Eq.(l)) at 
density (n) = 1/2 obtained from the analysis of a chain with 16 sites. The region 
kp > 1 is where superconductivity dominates in the ground state. When kp < 1 and 
at small Vcoui, SDW and CDW competes . At large values of Vcoui, or to the right of 
the superconducting region, the ground state is dominated by CDW order (states I, 
II, III, ... having different charge modulations). The open (full) circles denote holes 
(electrons). The size of each microscopic cluster in the CDW state increases as J/t 
increases, and the many CDW phases after IV are not shown. The dashed lines are 
tentative boundaries separating the different charge orders. Phase separation (PS) is 
only stable at Vcoui = in the bulk limit. 
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3. (a) kp in the ID t-J model, at density (n) = 1/2 on a 16 sites chain. The sohd 
hne corresponds to Vcoui = 0.00 and shows the divergence of kp at the point of phase 
separation. On the other hand at Vcoui = 0.15t, kp remains finite; (b) Behavior of 
the pairing and charge susceptibihties (xsup and Xcdw, respectively) vs J/t at fixed 
Vcoui = O.lOt; (c) Phase diagram of the 2D t-J model at density (n) = 1/2 obtained 
using a 16 sites square cluster. SC is (d-wave) superconductivity, PM a paramagnetic 
region, PS phase separation, I is the standard Wigner crystal, II a Wigner crystal of 
pairs, and III and IV are shown in the figure. Note the strong similarities with Fig. 2. 

4. (a) Mean-field phase diagram of the 2D t-U-V model including 1/r interactions at half- 
filling. PS denotes phase separation (stable over a finite region since in this calculation 
since the 1/r was truncated at four lattice spacings), d — wave is the superconducting 
phase, while FL denotes a Fermi-liquid region obtained by turning off the supercon- 
ducting order parameter in the mean-field calculation. I is a CDW Wigner crystal, II 
is a one column striped phase (i.e. phase III of Fig. 3c), and III is shown; (b) Generic 
representation of the phase diagram of electronic models with competing short-range 
attractive and long-range repulsive forces. For details see the text. 
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